Grobstruktursimulationen zur Untersuchung des Einflusses von Landoberflächenheterogenitäten auf die Struktur der konvektiven Grenzschicht sowie Analyse der Repräsentativität von Turbulenzmessungen

[no abstract

Heterogeneity-Induced Heat-Flux Patterns in the Convective Boundary Layer: Can they be Detected from Observations and is There a Blending Height? - A Large-Eddy Simulation Study for the LITFASS-2003 Experiment

An understanding of how the convective boundary layer (CBL) is mixed under heterogeneous surface forcing is crucial for the interpretation of area-averaged turbulence measurements. To determine the height and degree to which a complex heterogeneous surface affects the CBL, large-eddy simulations (LES) for two days of the LITFASS-2003 experiment representing two different wind regimes were undertaken. Spatially-lagged correlation analysis revealed the turbulent heat fluxes to be dependent on the prescribed surface flux pattern throughout the entire CBL including the entrainment layer. These findings prompted the question of whether signals induced by surface heterogeneity can be measured by airborne systems. To examine this question, an ensemble of virtual flights was conducted using LES, according to Helipod flight measurements made during LITFASS-2003. The resulting ensemble-averaged heat fluxes indicated a clear dependence on the underlying surface up to the top of the CBL. However, a large scatter between the flux measurements in different ensemble runs was observed, which was the result of insufficient sampling of the largest turbulent eddies. The random and systematic errors based on the integral length scale did not indicate such a large scatter. For the given flight leg lengths, at least 10-15 statistically independent flight measurements were necessary to give a significant estimate of heterogeneity-induced signals in the CBL. The need for ensemble averaging suggests that the observed blending of heterogeneity-induced signals in the CBL can be partly attributed to insufficient averaging.DFG/RA 617/21-

How Does the Choice of the Lower Boundary Conditions in Large-Eddy Simulations Affect the Development of Dispersive Fluxes Near the Surface?

Large-eddy simulations (LES) are an important tool for investigating the longstanding energy-balance-closure problem, as they provide continuous, spatially-distributed information about turbulent flow at a high temporal resolution. Former LES studies reproduced an energy-balance gap similar to the observations in the field typically amounting to 10–30% for heights on the order of 100 m in convective boundary layers even above homogeneous surfaces. The underestimation is caused by dispersive fluxes associated with large-scale turbulent organized structures that are not captured by single-tower measurements. However, the gap typically vanishes near the surface, i.e. at typical eddy-covariance measurement heights below 20 m, contrary to the findings from field measurements. In this study, we aim to find a LES set-up that can represent the correct magnitude of the energy-balance gap close to the surface. Therefore, we use a nested two-way coupled LES, with a fine grid that allows us to resolve fluxes and atmospheric structures at typical eddy-covariance measurement heights of 20 m. Under different stability regimes we compare three different options for lower boundary conditions featuring grassland and forest surfaces, i.e. (1) prescribed surface fluxes, (2) a land-surface model, and (3) a land-surface model in combination with a resolved canopy. We show that the use of prescribed surface fluxes and a land-surface model yields similar dispersive heat fluxes that are very small near the vegetation top for both grassland and forest surfaces. However, with the resolved forest canopy, dispersive heat fluxes are clearly larger, which we explain by a clear impact of the resolved canopy on the relationship between variance and flux–variance similarity functions. © 2021, The Author(s)

Modeling of land–surface interactions in the PALM model system 6.0: land surface model description, first evaluation, and sensitivity to model parameters

In this paper the land surface model embedded in the PALM model system is described and evaluated against in situ measurements at Cabauw, Netherlands. A total of 2 consecutive clear-sky days are simulated, and the components of surface energy balance, potential temperature, humidity, and horizontal wind speed are compared to observations. For the simulated period, components of the energy balance are consistent with daytime and nighttime observations, and the daytime Bowen ratio also agrees fairly well with observations. The model simulates a more stably stratified nocturnal boundary layer than the observations, and near-surface potential temperature and humidity agree fairly well during the day. Moreover, a sensitivity analysis is performed to investigate dependence of the model on land surface and soil specifications, as well as atmospheric initial conditions, because they represent a major source of uncertainty in the simulation setup. It is found that an inaccurate estimation of leaf area index, albedo, or initial humidity causes a significant misrepresentation of the daytime turbulent sensible and latent heat fluxes. During the night, the boundary-layer characteristics are primarily affected by surface roughness and the applied radiation schemes.publishedVersio

On the Effect of Surface Heat-Flux Heterogeneities on the Mixed-Layer-Top Entrainment

We used a set of large-eddy simulations to investigate the effect of one-dimensional stripe-like surface heat-flux heterogeneities on mixed-layer top entrainment. The profiles of sensible heat flux and the temporal evolution of the boundary-layer depth revealed decreased entrainment for small heat-flux amplitudes and increased entrainment for large heat-flux amplitudes, compared to the homogeneously-heated mixed layer. For large heat-flux amplitudes the largest entrainment was observed for patch sizes in the order of the boundary-layer depth, while for significantly smaller or larger patch sizes entrainment was similar as in the homogeneous case. In order to understand the underlying physics of this impact, a new approach was developed to infer local information on entrainment by means of the local flux divergence. We found an entrainment maximum over the centre of the stronger heated surface patch, where thermal energy is accumulated by the secondary circulation (SC) that was induced by the surface heterogeneity. Furthermore, we observed an entrainment maximum over the less heated patch as well, which we suppose is to be linked to the SC-induced horizontal flow convergence at the top of the convective boundary layer (CBL). For small heat-flux amplitudes a counteracting effect dominates that decreases entrainment, which we suppose is the horizontal advection of cold air in the lower, and warm air in the upper, CBL by the SC, stabilizing the CBL and thus weakening thermal convection. Moreover, we found that a mean wind can reduce the heterogeneity-induced impact on entrainment. If the flow is aligned perpendicular to the border between the differentially-heated patches, the SC and thus its impact on entrainment vanishes due to increased horizontal mixing, even for moderate wind speeds. However, if the flow is directed parallel to the border between the differentially-heated patches, the SC and thus its impact on entrainment persists.DFG/RA 617/21-1DFG/RA 617/20-1Niedersächsische Technische Hochschule (NTH

Importance of radiative transfer processes in urban climate models: A study based on the PALM 6.0 model system

Including radiative transfer processes within the urban canopy layer into microscale urban climate models (UCMs) is essential to obtain realistic model results. These processes include the interaction of buildings and vegetation with shortwave and longwave radiation, thermal emission, and radiation reflections. They contribute differently to the radiation budget of urban surfaces. Each process requires different computational resources and physical data for the urban elements. This study investigates how much detail modellers should include to parameterize radiative transfer in microscale building-resolving UCMs. To that end, we introduce a stepwise parameterization method to the Parallelized Large-eddy Simulation Model (PALM) system 6.0 to quantify individually the effects of the main radiative transfer processes on the radiation budget and on the flow field. We quantify numerical simulations of both simple and realistic urban configurations to identify the major and the minor effects of radiative transfer processes on the radiation budget. The study shows that processes such as surface and vegetation interaction with shortwave and longwave radiation will have major effects, while a process such as multiple reflections will have minor effects. The study also shows that radiative transfer processes within the canopy layer implicitly affect the incoming radiation since the radiative transfer model is coupled to the radiation model. The flow field changes considerably in response to the radiative transfer processes included in the model. The study identified those processes which are essentially needed to assure acceptable quality of the flow field. These processes are receiving radiation from atmosphere based on the sky-view factors, interaction of urban vegetation with radiation, radiative transfer among urban surfaces, and considering at least single reflection of radiation. Omitting any of these processes may lead to high uncertainties in the model results.publishedVersio

High-resolution large-eddy simulation of indoor turbulence and its effect on airborne transmission of respiratory pathogens - Model validation and infection probability analysis

High-resolution large-eddy simulation (LES) is exploited to study indoor air turbulence and its effect on the dispersion of respiratory virus-laden aerosols and subsequent transmission risks. The LES modeling is carried out with unprecedented accuracy and subsequent analysis with novel mathematical robustness. To substantiate the physical relevance of the LES model under realistic ventilation conditions, a set of experimental aerosol concentration measurements are carried out, and their results are used to successfully validate the LES model results. The obtained LES dispersion results are subjected to pathogen exposure and infection probability analysis in accordance with the Wells-Riley model, which is here mathematically extended to rely on LES-based space- and time-dependent concentration fields. The methodology is applied to assess two dissimilar approaches to reduce transmission risks: a strategy to augment the indoor ventilation capacity with portable air purifiers and a strategy to utilize partitioning by exploiting portable space dividers. The LES results show that use of air purifiers leads to greater reduction in absolute risks compared to the analytical Wells-Riley model, which fails to predict the original risk level. However, the two models do agree on the relative risk reduction. The spatial partitioning strategy is demonstrated to have an undesirable effect when employed without other measures, but may yield desirable outcomes with targeted air purifier units. The study highlights the importance of employing accurate indoor turbulence modeling when evaluating different risk-reduction strategies

The Parallelized Large-Eddy Simulation Model (PALM) version 4.0 for atmospheric and oceanic flows: model formulation, recent developments, and future perspectives

In this paper we present the current version of the Parallelized Large-Eddy Simulation Model (PALM) whose core has been developed at the Institute of Meteorology and Climatology at Leibniz Universität Hannover (Germany). PALM is a Fortran 95-based 5 code with some Fortran 2003 extensions and has been applied for the simulation of a variety of atmospheric and oceanic boundary layers for more than 15 years. PALM is optimized for use on massively parallel computer architectures and was recently ported to general-purpose graphics processing units. In the present paper we give a detailed description of the current version of the model and its features, such as an embedded 10 Lagrangian cloud model and the possibility to use Cartesian topography. Moreover, we discuss recent model developments and future perspectives for LES applications.DFG/RA/617/3DFG/RA/617/6DFG/RA/617/16DFG/RA/617/27-

Development of an atmospheric chemistry model coupled to the PALM model system 6.0: Implementation and first applications

In this article we describe the implementation of an online-coupled gas-phase chemistry model in the turbulence-resolving PALM model system 6.0 (formerly an abbreviation for Parallelized Large-eddy Simulation Model and now an independent name). The new chemistry model is implemented in the PALM model as part of the PALM-4U (PALM for urban applications) components, which are designed for application of the PALM model in the urban environment (Maronga et al., 2020). The latest version of the Kinetic PreProcessor (KPP, 2.2.3) has been utilized for the numerical integration of gas-phase chemical reactions. A number of tropospheric gas-phase chemistry mechanisms of different complexity have been implemented ranging from the photostationary state (PHSTAT) to mechanisms with a strongly simplified volatile organic compound (VOC) chemistry (e.g. the SMOG mechanism from KPP) and the Carbon Bond Mechanism 4 (CBM4; Gery et al., 1989), which includes a more comprehensive, but still simplified VOC chemistry. Further mechanisms can also be easily added by the user. In this work, we provide a detailed description of the chemistry model, its structure and input requirements along with its various features and limitations. A case study is presented to demonstrate the application of the new chemistry model in the urban environment. The computation domain of the case study comprises part of Berlin, Germany. Emissions are considered using street-type-dependent emission factors from traffic sources. Three chemical mechanisms of varying complexity and one no-reaction (passive) case have been applied, and results are compared with observations from two permanent air quality stations in Berlin that fall within the computation domain. Even though the feedback of the model's aerosol concentrations on meteorology is not yet considered in the current version of the model, the results show the importance of online photochemistry and dispersion of air pollutants in the urban boundary layer for high spatial and temporal resolutions. The simulated NOx and O3 species show reasonable agreement with observations. The agreement is better during midday and poorest during the evening transition hours and at night. The CBM4 and SMOG mechanisms show better agreement with observations than the steady-state PHSTAT mechanism. © 2021 Copernicus GmbH. All rights reserved

Sensitivity analysis of the PALM model system 6.0 in the urban environment

Sensitivity of the PALM model 6.0 with respect to land-surface and building properties is tested in a real urban environment in the vicinity of a typical crossroads in a densely built-up residential area in Prague, Czech Republic. The turbulence-resolving PALM is able to simulate the urban boundary layer flow for realistic setups. Besides an accurate representation of the relevant physical processes, the model performance also depends on the input data describing the urban setup, namely the building and land-surface properties. Two types of scenario are employed. The first one is the synthetic scenarios altering mainly surface and material parameters such as albedo, emissivity or wall conductivity, testing sensitivity of the model simulations to potentially erroneous input data. Second, urbanistic-type scenarios are analysed, in which commonly considered urban heat island mitigation measures such as greening of the streets or changing surface materials are applied in order to assess the limits of the effects of a particular type of scenario. For the synthetic scenarios, surface parameters used in radiation balance equations are found to be the most sensitive overall followed by the volumetric heat capacity and thermal conductivity of walls. Other parameters show a limited average effect; however, some can still be significant during some parts of the day, such as surface roughness in the morning hours. The second type, the urbanistic scenarios, shows urban vegetation to be the most effective measure, especially when considering both physical and biophysical temperature indicators. The influence of both types of scenario was also tested for air quality, specifically PM2.5 dispersion, which generally shows opposite behaviour to that of thermal indicators; i.e. improved thermal comfort brings deterioration of PM2.5 concentrations. © 2021 Michal Belda et al